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Abstract. Measurements by Ulysses have confirmed 
that there is no significant gradient with respect to he- 
liomagnet ic latitude in the radial component, B r , of the 
interplanetary magnetic field. In the corona, the plasma 
{3 is << 1, except directly above streamers, so longitudi- 
nal and latitudinal gradients in field strength will relax 
due to the transverse magnetic pressure gradient force 
as the solar wind carries magnetic flux away from the 
Sun. This happens quickly enough so that the field is 
essentially uniform by 5 — 10 Rq, apparently remaining 
so as it is carried to beyond 1 AU. Here, we illustrate the 
coronal relaxation with a qualitative physical argument 
and by reference to a detailed MHD simulation. 

1. Introduction 

Ulysses in 1993-1995 [Smith and Balogh, 1995; Balogh 
et al, 1995] and ICE and IMP-8 in 1984-1988 [Burton 
et al , 1995] observed no significant gradient in helio- 
magnetic latitude in the radial component, B r , of the 
interplanetary magnetic field (IMF) at heliocentric dis- 
tances of 1-4 AU. These observations were made near 
solar minimum when the Sun’s magnetic field is most 
nearly an axially aligned dipole. Figure 1 shows the 
results from IMP-8 and ICE scaled to 1 AU and plot- 
ted against magnetic latitude. The data are five degree 
bin averages. Curves (a), (b), and (c) are from mod- 
els which will be discussed below. The apparent small 
gradient in B r near the magnetic equator is probably 
due to small errors in sector identification so that the 
true gradient is completely negligible. Essentially no 
gradient was observed by Ulysses in 1993-1995 up to a 
heliographic latitude of 80°. Here we discuss the rea- 
son why the magnetic field at the top of the corona, 
i.e. at 10 Rq, does not have any significant gradient in 
latitude or longitude outside of the heliospheric plasma 
sheet (HPS) [Gosling et al, 1981] surrounding the he- 
liospheric current sheet (HCS). In §2, we explain why 
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the coronal imprint on the IMF lacks latitudinal and 
longitudinal (i.e. (0, <f>)) gradients on a physical basis. 
This is verified in Section 3 using results from an MHD 
model of the corona. In a second study, we examine 
the absence of further redistribution of magnetic flux in 
the interplanetary medium by latitudinal gradients in 
thermal pressure [Suess et al, 1996]. 

2. The Coronal Magnetic Field 

The immediate reason for the absence of ( 0 , $ ) gradi- 
ents in the magnetic field at 10 Rq is that the plasma (3 
(the thermal pressure divided by the magnetic pressure) 
is small between 1.0 Rq and ~ 10.0 Rq except in the 
HPS above streamers. This is demonstrated here using 
a potential field-source surface (PFSS) model of coronal 
structure as a tool. PFSS models are routinely used to 
estimate coronal streamer locations and sizes and the 
location of the magnetic neutral line lying above the 
streamer belt, which defines the base of the heliospheric 
current sheet (HCS) [Hoeksema, 1991]. However, they 
fail to predict the absence of ( 0 , <j>) gradients in the in- 
terplanetary medium. In spite of this failure, a PFSS 
model can show that the gradients are smoothed out 
over a very short distance inside 10 Rq. This is done 
by using the PFSS model to demonstrate its own inter- 
nal inconsistency and to estimate the speed at which 
the gradients are therefore smoothed. 

In a PFSS model the field is assumed to be current 
free between the photosphere and a “source surface” 
placed at about 2.5 Rq , the source surface is assumed 
spherical, and the magnetic field is defined to be ra- 
dial at the source surface. The boundary condition at 
1 .0 Rq is that a scalar be specified over the sphere - nor- 
mally the observed line-of-sight field. Figure 2 shows 
the line of sight photospheric field and the resulting 
source surface magnetic field in a PFSS model using 
Wilcox Solar Observatory data for Carrington Rotation 
1779 (August-September 1986) [Hoeksema, 1991]. This 
rotation is included in the time period of IMP-8 and 
ICE data used in Figure 1, and it is obvious the source 
surface field varies smoothly up to and across the neu- 
tral line, unlike the IMF which exhibits a discontinuity 
at the HCS. The implied transverse gradients in | B 2 |, 
d | B 2 | fd<j> and d | B 2 | jdO , which are completely 
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Figure 1 . Radial magnetic field values in five degree 
latitude bins, scaled to 1 AU, from ICE and IMP-8 dur- 
ing 1984-1988, bracketing solar minimum [Burton et ai 
1995]. (a) is a pure current sheet model [Wolfeon, 1985], 
at a heliocentric distance of more than 50 Rq. (b) is 
from Pneuman and Kopp’s [1971] MHD model, at a 
large distance from the base, (c) is the field variation 
for a potential field source surface model with the source 
surface at 2.5 Rq. 


inconsistent with the observations shown in Figure 1, 
would have to be balanced by corresponding gradients 
in the pressure to persist. 

The distance it will take for the gradients to relax 
can be calculated using the PFSS model to estimate 
the relaxation speed. Assume the coronal magnetic field 
between 1 Rq and the source surface at R s to be defined 

by 

B = -V$ V 2 $ = 0 (1) 

Taking the field to be radial on R 3 , axisymmetric, 
and dipolar at 1 Rq where B r (R Q ) = B ro cosQ — 
B ro P® (cos 0) gives 
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A plot of the field lines from (2) is shown in Figure 
3. The solid lines are field lines derived from (2) and 
the dashed lines are for a vacuum dipole in which all 
of the field lines return to the base. For the physical 
parameters, we use the observed radial field strength at 
1 AU of 3.5 x 10 -5 G [Smith and Balogh, 1995] and 
extrapolate it back to the Sun to give B ro — 4.4 G. 
For the temperature and density, let T(R S ) — 10 6 K 
and n(R s ) = 3 x 10 5 cm -3 [Allen, 1964; Habbal et ai , 
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Figure 2. Synoptic charts of the photospheric (bottom) and computed coronal (top) magnetic 
fields in September 1986 - Carrington Rotation 1779. The heavy solid fine shows the neutral 
line dividing inward and outward polarity regions. The dashed contours and heavy shading show 
negative (inward) field regions. Solid contours and light shading show positive areas. Shading 
begins at ±200 //T at the source surface located at 2.5 Rq . No data above 70° are shown 
because they are not measured independently due to the low spatial resolution of the Wilcox 
Solar Observatory observations [Hoeksema, 1991]. 







Figure 3. Magnetic field lines for a vacuum dipole 
(dashed lines) and a PFSS model with the source sur- 
face placed at 2.5 Rq (solid lines) in which the bound- 
ary condition on the source surface is that the field be 
radial. The source surface is at the outer half-circle. 


1995]. Then, with (3 = 2nfeT/(J3^/87r), evaluation on 
the source surface gives 

B r (R s ,0) = 0.4 cos 0 G (3) 

13(6, R s ) = 0.0126 cos -2 0 (4) 

such that (3(6 — 0, R s ) — 0.0126 and j3(0 = 80°, R a ) = 
0.42. Curve (c) in Figure 1 shows equation (2) extrap- 
olated to 1 AU and compared with the observed field, 
again illustrating the failure of the PFSS model to pre- 
dict the observed variation. Equation (3) shows that 
for typical parameters, f3 « 1 over the bulk of the 
open region. Only directly over the streamer is (3 > 1, 
asymptotically approaching oo as 9 — ► 90°. Inside the 
source surface, (3 decreases further, while it grows only 
slowly until it becomes of 0[1] only in the vicinity of 
10 Rq. The implied transverse magnetic pressure im- 
balance cannot be sustained by thermal pressure gradi- 
ents and will therefore be reduced, or relax, in the solar 
wind frame of reference. This occurs approximately at 
the fast mode speed, v/ 7 because the magnetic pressure 
gradients and the consequent relaxation are perpendic- 
ular to the direction of the magnetic field. 

The fast mode speed is given by 


, 2 kT B*(R.,6) 

V f {()) = ^ + 47mm p 


(5) 


such that Vf(0) « 1620 cos 0 km/s. Assume for this 
example that the radial flow speed of the solar wind in 
the mid-corona is v r — 400 km/s. This is probably large 
compared to typical coronal flow speeds and therefore 
will give an upper bound to the relaxation distance, 
A r . For a tranverse relaxation through a latitude angle 
(6 — 6 3 ) , at speed vj , in solar wind moving at speed v r , 


\ r can be estimated by using the angular speed given 
by 0 — (vf/r) = dO/dt. Substituting v r — dr/dt and 
integrating gives 

d9 = ( v r /r)dt — (vj fr)(dr/v r ) (6) 

0 - 9 S = (vf/v r ) In (r/R a ) (7) 

where, A r = (r — R s ). 0 S is the starting angle at R s . 
Taking (0 — 0 3 ) = 45° = 7r/4 radians to represent a 
large transverse angular displacement and assuming a 
nominal mean value for Vf of 1200 km/s, then A r — 
0.3 R s = 0.75 Rq. The relaxation is therefore complete 
inside 5 Rq. This could have been anticipated because 
v/ ^ v r in this region. So, rather than (2) representing 
the coronal imprint on the IMF radial magnetic field 
strength, it is appropriate to take the field strength as 
constant at the top of the corona everywhere except in 
the HPS. 

3. Comparison With An MHD Model 

The above argument based on the fast mode wave 
speed in the corona is fully supported by self-consistent 
MHD calculations of the relaxation of the transverse, 
gradients [e.g., Suess, et ai 1977, Fig. 2; Steinolfson, 
Suess, and Wu, 1982; Steinolfson, 1988]. The results 
from the earliest such calculation [Pneuman and Kopp, 
1971] are shown by curve (b) in Figure 1, where there 
is a large, but not complete relaxation in the transverse 
field gradients - probably being incomplete due to the 
photospheric field strength in this calculation being only 
1 G. A full relaxation is explicitly shown in the quan- 
titative global coronal model of Wang, Wu, Suess, and 
Poletto [1995] which we use here as an example. The 
magnetic field variation with polar angle for the steady 
state is shown in Figure 4. This model has a steady 
volumetric heat source that varies with radius and po- 
lar angle, being smaller inside the equatorial streamer 
than in the polar coronal hole, and a field strength at 
the base of 0[3] G. It produces densities in the hole 
and streamer of the same magnitudes as observed. The 
location of the streamer and the HPS lying above the 
streamer, where (3 > 1, are shown as the shaded re- 
gion in Figure 4. This figure explicitly shows that the 
field strength is already nearly constant by 2.47 Rq , 
and by 4.88 Rq the variation is insignificant - outside 
the HPS. This result is a detailed confirmation of the 
estimates made in §2 and closely reproduces the obser- 
vations shown in Figure 1. The model further shows 
that (3 « 1 everywhere inside ~ 10 Rq except in the 
HPS. Beyond ~ 10 Rq , (3 > 1. 

Returning to PFSS models, these qualitatively simu- 
late the field strength variation low in the corona, inside 
~ 2.0 Rq. Conversely, a simple current sheet model of 
the corona leads to a uniform field strength far from 
the Sun as shown by curve (b) in Figure 1 [Wolfson, 
1985]. Combining external current sheets with PFSS 
models has been successful in simulating the dynamics 
of the relaxation of the transverse magnetic field gra- 
dients [Zhao and Hoeksema, 1995; Wang and Sheeley, 
1995] and is a way of retaining the practical applications 
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Figure 4. The total magnetic field strength at: (a) 

1.00 i?0, (b) 1.47 Rq, (c) 2.47 Rq , (d) 4.88 Rq versus 
polar angle. These values are taken from the global 
coronal simulation of Wang, Wu, Suess, and Poletto 
[1995] in which there is a volumetric heat source that 
produces a realistically thin heliospheric plasma sheet. 
The solid lines are the model described here. 

of these relatively simple models. However, it must be 
recognized that they do not do a good job of modeling 
the transverse variation of the field strength between 

2.0 and 10.0 Rq. 

Finally, the Ulysses, ICE, and IMP-8 results imply 
that latitudinal thermal pressure gradients are insuffi- 
cient to redistribute magnetic flux in the interplanetary 
medium between 10 Rq and 1-4 AU. We have analyzed 
the latitudinal temperature and density gradients ob- 
served at Ulysses and combined this with an analytic 
calculation of the resulting flux redistribution to show 
that an approximately constant field strength at 1 AU 
is consistent with the observed pressure gradients [Suess 
et ai, 1996]. 
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